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Abstract. 

Nonperturbative dynamic theory has a particular advantage in studying the 
transport in a quantum impurity system in a steady state. Here, we develop a new 
approach for obtaining the retarded Green's function expressed in resolvent form. We 
use the Heisenberg picture to facilitate dynamic theory and propose a new systematic 
method of collecting the basis vectors spanning the Liouville space, which is the most 
crucial step in obtaining the resolvent Green's function. We obtain all the linearly 
independent basis vectors for studying the single-impurity Anderson models with 
one and two reservoirs. The latter is an appropriate model for studying the Kondo 
phenomenon in a steady state when a bias is applied. This is one of long standing 
subjects in theoretical condensed matter physics. 
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The single-particle Green's function is a basic tool for studying the correlation 
effects in many-body systems. However, calculation of this function for a strongly 
correlated system is not usually successful because nonperturbative treatment is 
generally required. Moreover, the situation becomes even worse for strongly correlated 
systems under a steady state. Such systems are quantum point contact [H El El H] , single 
electron transistor with a quantum dot J5j El [TJ [8] , and scanning tunneling microscopy 
on a magnetic atom adsorbed on a metallic surface [91 [101 [Til EG] • These steady state 
Kondo phenomena remain unexplored because no successful theory, which can handle 
the nonequilibrium steady state for a quantum impurity system with strong electron 
correlation, is available. 

A well-known formula for a steady current through a small interacting region 
connected to charge reservoirs having different chemical potentials has been presented 
more than a decade ago [TBI E] . The formula is given by 

e rdu T L (u)T R (u) . 
J=-J-J ^ TL{uj) + TR{u) ihH - ImG+» (1) 

for proportional lead functions, i.e., T L (oo) oc T r (oj), where L and R denote the left and 
right reservoirs, respectively and f(oo) denotes the Fermi distribution. The retarded 
Green's function G^ da (uj) for the mediating atom must be obtained for a given bias. 
However, this steady state Green's function has never been obtained for an atom at 
which a strong interaction exists. 

The purpose of this study is to present a new method of determining the steady 
state Green's function appearing in Eq. (JJJ). For this purpose, we adopt the resolvent 
Green's function in the Heisenberg picture |15j . which is given by for a fermion system 
as 

iGt Ja {z) = (c la \{zl + il,)- 1 \c 3a ), (2) 

where z = —iuj + r\ and q CT indicates a fermion of spin o annihilating at a state or site %. 
The symbols I and L in Eq. (J2J) denote the identity operator and the Liouville operator, 
respectively. The latter is defined by LO = [H, O], where TC and O are the Hamiltonian 
and an arbitrary operator, respectively. Although we are familiar with the resolvent 
form in the Schrodinger picture, 

G+(w) = <&|(w + ir ] -H)- 1 \cf> j ), (3) 

on the basis of the static basis {<f>j\j = 0, 1, - • • , oo}, we employ Eq. (T5]) using 
the operator basis in this study because dynamical approach is more appropriate for 
handling steady state situations. 

Calculating the retarded Green's function in resolvent form starts from ensuring a 
complete set of basis vectors. A notable point is that the resolvent form does not provide 
information about the basis vectors. Therefore, the Lanczos algorithm [161 E! m the 
Schodinger picture and the projection operator technique [TU [19] in the Heisenberg 
picture are used to obtain the static and dynamic basis vectors, respectively. These 
techniques, however, require a heavy work for a nontrivial system because they employ 
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the Gram-Schmidt orthogonalization procedure for the linearly independent vectors 
7i n |0o) or L n Cjo-, where n — 0, 1, • • • , oo. For this reason, we propose a new methodology 
that provides a simpler procedure to obtain the basis vectors. 

The retarded Green's function G dd(T (uj) can be obtained by calculating the matrix 
inverse (M -1 )^, where M is composed of the elements 

My = zSij - (iLej\e,i) = z5ij + (e^Le^), (4) 

once a complete basis {£i\£ = l,---,oo} spanning the Liouville space is given. We 
choose e± = c^, to obtain G dda (u). The inner product in Eq. (0J is defined by 
(efc|eg) = ({ek,e\}), where the angular and curly brackets denote statistical average 
and anticommutator, respectively. We now present a systematic method to collect the 
basis vectors. 

In the previous report [20J, we have obtained the basis vectors intuitively. In 
this study, however, we provide a systematic manner to secure universal validity and 
extensive applicability. Expressing the resolvent Green's function in a matrix form such 

as 

*s,w - «4 m ( I ) ( [J; ) (5) 

is the most crucial step to obtain the basis vectors in systematic manner. Once we get 
this form for the Green's function, it is convinced that the Liouville space of Gfj a (z) 
will be spanned completely by the linearly independent components of vectors \A) and 
\B). 

To obtain the matrix product form in Eq. (jSJ), we expand the Green's function 
operator G = [(zl + iLj) + iLc}^ 1 , where Lj and Lc represent the Liouville operator 
for the isolated part of the Hamiltonian Tii and the connecting part Tic, respectivel, in 
powers of Lc using the operator identity (A + B)^ 1 = A^ 1 — A~ l B(A + B)^ 1 . Then, 
one can express the retarded Green's function as follows: 

iGf ja (z) = (c i0 .|G z |cj>) - (cirlGiliLcGiCjr) + {c ia G I iL c \G\i'LcG I Cj <T ), (6) 

where Gi = (zl + iLi)~ l . Equation ([6l) is transformed into a matrix form, 

iG+ a (z) = ({c i(T \ ($i|)G(|c^> |$,)) T , (7) 

(G — G \ -> 
^ ~ j, = \ihcGiCj a ), and the superscript T denotes the 

transpose. After a similarity transformation that diagonalizes G, Eq. ([7]) is given by 

iG+ a (z) = ((c i<T \ (^\)G D (\c ja ) |$,» T , (8) 

where G^ = ^ ^ q ^ and 

\cja) = \c jrT ) + (Gi - = \c jtT ) + \Lc\-izI + L)$j). (9) 

Equation (jHJ) is the form we want to have and the linearly independent components of 
vectors \cj^) and completely span the Liouville space of Gfj a (u). In conclusion, 
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the systematic method for collecting the basis vectors involves finding all linearly 
independent components comprising the vector \£j a ) because |<3> 3 -) is contained in it. 

We will subsequently demonstrate the determination of the basis vectors with an 
example. We are interested in a mesoscopic system with a mediating Kondo atom 
between two metallic reservoirs. This system can be described by a single-impurity 
Anderson model with two metallic reservoirs, whose Hamiltonian is written as 

H = Yl e k c kl c L + Yl e dC da c da + U n d ^n di 

k,a,a=L,R a 

+ e (Vkdcict + v k * d 4ic da ), 

k,a,a=L,R 

where e k , e d , V k d, and U indicate the energies of an electron of momentum k in the lead, 
the level of a mediating atom, the hybridization between the atom and lead, and the on- 
site Coulomb repulsion at the atom, respectively. Because the two-reservoir Anderson 
model is a straightforward extension of the one-reservoir model, we first consider a 
one-reservoir Anderson model. Its isolated and coupled parts are given by 

= £ € k C L C ka + e d C da°da + U U^H^ (10) 
k,a a 

and 

n c = £(^4 CfaT + v k * d clc da ). (ii) 

k,a 

We first demonstrate the calculation of Q d = iL c (J + £L/) -1 c<ff and then \c da ) for 
Hi and TCc above. One can easily see that (I + iL/) _1 Q| yields only two operators, c d ^ 
and n d ic d p Applying he to these operators changes the index d(k) into k(d) in fermion 
operators and n d i into j d ,. Therefore, the components of § d are given by 

$d = (cfc T , n dl c kh j^car), (12) 

where k — 1, 2, ■ • • , oo and j d y = i J2k(VkdC da c ka — V kd c k(T c d(T ). We now focus on the vector 
\c dr7 ) in Eq. The operator L$ d = [H,c k ^ + [H,n di c k ^} + [H,j di c d ^] yields c rfT and 
Cfcf from the first commutator, n d ic d ^, n d ic k ^, and j d ^c k ^ from the second commutator, 
and j dl c dh j dl c kh jJc dT , j dl n dl c dh and j dl n dr c d ^ from the third commutator. The 
last operator j d \n d -\c d -\ is a vanishing one, and the operator j d \n d ic d ^ is dynamically 
equivalent to jai c <t\- Thus, the linearly independent components of (I + L)$ rf are 
classified into two groups, i.e., one involving c k \ and the other involving c d p They 
are 

(I + L)$^ = (c feT , n dl c kh j dl c kr ), (13) 

for k — 1, 2, • • • , oo, and 

(I + L)$^ = (c dT , n dl c dh j dl c dh jJc dT ). (14) 

As a last step, we apply L^ 1 , which is equivalent to the repeated application of h c , 
to the operators in Eqs. ( |T3l) and ( |T4l) . One can easily see that the multiple application of 
he to Cfcf, Crff , and n d -\ simply reproduces the existing operators and additionally j d \C k ^. 
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Remaining linearly independent operators are come from Lg>j'J , where n = l,2,---,oo. 
Thus, one may classify the basis vectors into two groups such as 

L c 1 ( I + L )®d = {cfcT, n dl c kh j'^c^, jjjcfct, (L^jJ)c feT } (15) 

and 

L5 1 (I + L)^ = {q lJT , n di c dT , j^c dT , jJc dT , (LJjJ)^}. (16) 

The operators shown in Eqs. (ITBl and ({TBI describe all possible linearly independent 
manners of annihilating an up-spin at the site, d, of the mediating atom in time t. They 
completely span the Liouville space of c^(t) in other words. 

It is important to identify the meaning of the basis vectors (L^'J )cfe| and (L^j^c^ 
for later discussion. Application of L c to j% gives rise to 

Lcj7i = -*'E XX^kd c ii c ki + VkdV^c^cn) + |y M r4i c ^ 

l k k 

and 

LciJ = EE vw^* - 53 E ^4|C U . 

Ik lk 

These operators represent a round trip of a down-spin electron between the mediating 
atom and metallic lead. Therefore, depicts ra-time trip of a down-spin electron 

between the atom and lead without coming into contact with an up-spin electron at the 
atom. This hardly happens in reality. 

Even though all linearly independent members of a basis of the Liouville space 
for the single-reservoir Anderson model have been obtained, working with these full 
members yields unnecessary complication that cannot be easily overcome. One of the 
advantages of the present method is that the physical meaning of each basis vector is 
clear. Therefore, it is possible to eliminate some of basis vectors that are unimportant 
in describing the dynamics of the system at a particular parameter regime. A particular 
concern must be given to the basis vector n^c^. Since n 2 d \ = n d ^ this basis vector 
involves all higher orders of [/-processes. Therefore, one must eliminate n d ^c d ^ in 
studying Kondo regime unless all [/-processes are treated exactly. However, n^Qj 
will be the most important basis vector when we study a small U -regime. As mentioned 
above, multiple round trips without coming into contact with other spin electron at 
the mediating atom are practically rare; therefore, one may further eliminate the basis 
vectors (LcJdD^T- In other words, we assume that the down-spin electron annihilates 
once it arrives at the mediating atom. This assumption may be valid enough to describe 
the spin exchange interaction in the Kondo regime. 

Now, turn our attention to the basis vectors combined with Ck, i.e., those shown in 
Eq. (}T5|) . This part of the basis vectors is transformed into the self-energy under matrix 
reduction that is discussed later. We neglect the basis vectors (L^j^c^ by the same 
reason neglecting (L£v/J )c d ^ above. We also find that the role of the basis vectors j^yC^ 
is the same as that of n^ic^- However, their contribution to constructing self-energy is 
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negligible compared with that of n^c^- Therefore, we further eliminate j^c^f. Then, 
we finally obtain the reduce the Liouville space constructed by the basis vectors 

(cdh Jdi c dh Jdi c dl) and (c fcT , n dl c k] ), 

where k = 1, 2, • • • , oo for describing the dynamics of the single- impurity Anderson 
model in the Kondo regime. 

Since the resolvent form of G dd Jz) is given by the Laplace transform of the 
coefficient a 1 (t) in the expansion c d ^(t) = Y^i^if)^^ where ii = qj, one must have 
an equality ai(t) = G dd ^(t). This is valid when all the basis vectors except c#\ are 
orthogonal to c d -\. This orthogonality condition is satisfied by introducing the expression 
5 A = A — (A) for the operators n d -\ and j%. Thus, we finally obtain the reduced Liouville 
space spanned by the basis vectors 

{cdh 5 Jdi°dh ^jj[Cd\) and (c kh 5n dl c kr ), 

where k = 1, 2, • • • , oo, for the one-reservoir Anderson model in the Kondo regime. In 
practical calculations, we use the normalized forms of the basis vectors. The basis 
vectors given above are the same as those used in the previous study [20] in which we 
obtained them intuitively. 

The retarded Green's function G^(uj) is obtained by calculating the matrix 
inverse (M -1 )^ constructed by the elements of Eq. (J4j) in terms of the basis given 
above. Calculating the inverse of an infinite-dimensional matrix M is a nontrivial 
problem. For this purpose, we perform matrix reduction by using Lowdin's partitioning 
technique [211 122] . One can reduce the infinite dimensional matrix M to a 3 x 3 matrix 
for the basis given above. The reduced M 3x3 is given by 

Mj x3 = M* xS - M OOX 3(M^ oxoo )- 1 M 3 xoo, (17) 

where M 3x3 is the part constructed by the basis vectors (c^, 8jJ\C d ^, 5j d ^c d -]), M^ xoo 
by the basis vectors (c^, dn^c^), where k = 1,2, •••,oo, and and M 3xoo 

represent two corner parts of the matrix M. The second term of Eq. (TlTj) constructs the 
self-energy of the retarded Green's function given by iG dd Ju) = [(M 3x3 ) -1 ] n . We will 
discuss the specific results for the single impurity Anderson model in a separate paper. 

The Liouville space for the two-reservoir Anderson model is a straightforward 
extension of that of the single-reservoir model. The basis vectors are given by 

(Cfcf 5n dl c k^)i wnere k = l,2,---,oo, 

( s Jdi L °dh s Jdi Lc dh °dh Sjdfcdh s Jdi Rc di), and 
($ n diCk V c£ T ), where k = 1, 2, • • • , oo. 

The superscripts L and R denote the left and right metallic leads, respectively. This basis 
is used in explaining the experiments for the steady state Kondo phenomena occurring 
in quantum point contact or quantum wireJTJ^ElH], single electron transistor [3 [6j[7l [8], 
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and scanning tunneling microscopy [HI QUI Ell E2]- The reduced matrix in this case is a 
5x5 matrix, i.e., Mg x5 , which is given by 

M^ x5 = M^ x5 - M^ x5 (M^ X0O )- 1 M^ xoo - M£ x5 (CJ _1 ML. (18) 
where Mg x5 is constructed by the basis vectors 

and M^ xoo denotes the oo x oo matrix constructed by the basis vectors 
(c£i, Sridic^), where k = 1,2, ■■■,oo. The second and third terms of Eq. ({TBI) give 
rise to the self-energy of the retarded Green's function, iG dd Juj) = [(M.f xn )^ 1 ] 33 , of the 
two-reservoir Anderson model. 

This is the end of formalism for the retarded Green's function via the new approach. 
However, calculating the matrix elements is another laborious process to obtain the 
retarded Green's function in a steady state. We present those calculations in a separate 
work. Even though one calculate the matrix elements, a barrier for the retarded Green's 
function still remains because the elements of the matrix M^j contain the undetermined 
quantities (nd a ), (jda' R )i an d {jdo-' R )> which are determined by the retarded Green's 
function. They are expressed by G dda (id) as follows: 

/ \ 1 r /^) r » + /flMr» t r + i w rim 
M = -- r ^M + r«M ** G **M**> ( 19 ) 

1 r did 



(20) 



and 



(jt LiR) ) = r -f L{R) (cd)T L ^(u;)ReG+M, (21) 

where f (w) = T L (u)r R (u)/[T L (u) + T R (u)}. The expressions of Eqs. (pED and ([20]) are 
valid only when T L (u) oc T r (uj). Because of these formal relations, the retarded Green's 
function and the above quantities must be obtained by running a self-consistent loop 

<m> (0) , Ul L ' R ) (0) - < 0) M - - Kt> (0 \ <;?'*> (0) - 

in terms of Eqs. flU-flU)) and = [(Mg x6 )-V 

In sum, we have presented a systematic methodology for selecting basis vectors 
spanning the Liouville space, which is the most crucial step in calculating the retarded 
Green's function in resolvent form, and have suggested a procedure for calculating the 
retarded Green's function. The method presented in this study has several advantages: 
(i) one can find a complete set of basis vectors systematically; (ii) the physical meanings 
of the basis vectors obtained by this method are very apparent, allowing one to identify 
and remove the unimportant basis vectors for a particular parameter regime and 
construct a reduced Liouville space; and (iii) the self-consistent calculation does not 
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require a nonequilibrium density matrix, which makes nonequilibrium problem difficult. 
As a final comment, we mention on the remaining unknowns in addition to (n^), (jd^' R ), 
and {jda ,R ) m the matrix elements. These unknowns may be used as free parameters 
characterizing the system under consideration. In a separate study, we present the 
spectral function and differential conductance for a specific system showing steady state 
Kondo phenomenon. 
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